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Abstract. We propose a method to calculate the two-dimensional (2D) fission-fragment yield Y[Z, N) 
versus both proton and neutron number, with inclusion of odd-even staggering effects in both variables. 

The approach is to use Brownian shape-motion on a macroscopic-microscopic potential-energy surface 
which, for a particular compound system is calculated versus four shape variables: elongation (quadrupole 
moment Q 2 ), neck d, left nascent fragment spheroidal deformation cfi, right nascent fragment deformation 
ef 2 and two asymmetry variables, namely proton and neutron numbers in each of the two fragments. The 
extension of previous models 1) introduces a method to calculate this generalized potential-energy function 
and 2) allows the correlated transfer of nucleon pairs in one step, in addition to sequential transfer. In the 
previous version the potential energy was calculated as a function of Z and N of the compound system 
and its shape, including the asymmetry of the shape. We outline here how to generalize the model from 
the “compound-system” model to a model where the emerging fragment proton and neutron numbers also 
enter, over and above the compound system composition. 

PACS. 25.85.Ca discribing text of that key - 24.10.Lx discribing text of that key - 24.75.-|-i discribing 
text of that key - 25.85.Jg discribing text of that key 


1 Introduction 

In previous investigations it has been shown that a re¬ 
alistic description of the experimentally observed fission- 
fragment charge distributions can be obtained by means of 
random walks on tabulated five-dimensional (5D) potential- 
energy surfaces calculated for a densely spaced grid for 
over five million different shapes [T1[^I5] . It was partic¬ 
ularly encouraging that the 70 charge-yield distributions 
measured at GSI [4] were well reproduced, including the 
transition from symmetric fission for light Th isotopes to 
asymmetric fission for the heavier isotopes beyond A ss 
226. However, the often strongly fluctuating patterns of 
odd-even staggering were not possible to obtain in that 
implementation of the model. Briefly stated, the reason 
was that there was no mechanism that allowed the prop¬ 
erties of the two nascent fragments to affect the calculated 
potential-energy surface. For example, pairing was treated 
for the compound system as a whole. 

In the literature there have been numerous discussions 
of odd-even staggering and models proposed to quanti- 
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tatively describe this feature. For example the effect has 
been correlated with Coulomb-related quantities, such as 

jA and Z'^ jA^/^, referred to as order parameters, and 
to pairing effects on the nuclear level density and its vari¬ 
ation with excitation energy, see Refs. [5116117] and refer¬ 
ences therein. Other models are based on properties of 
separated fragments and thermal equilibrium at scission. 
For a review see Ref. [5]. Common for these models are 
that they are not based on detailed, calculated potential- 
energy surfaces or dynamical evolution on such surfaces. 
The models often also contain a substantial number of pos¬ 
tulated terms with parameters which are determined from 
adjustments to observed yields. Another group of models 
do treat dynamical evolution, in a Langevin approach. Un¬ 
til now they are based on macroscopic potential-energy 
surfaces, or, when shell effects are included the calcula¬ 
tions are performed for three shape variables for fairly 
high excitation energies {E* = 20 MeV), see Ref. |5] for a 
brief review but with extensive references to original work. 

In the Brownian shape-motion (BSM) model the only 
parameters are those of a well-established macroscopic- 
microscopic model used to calculate the five-dimensional 
(5D) potential-energy surfaces |1()] in which the parame- 
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ters have been unchanged since 2002 the additional 
critical neck radius at which the mass split is frozen, a 
weak bias potential (but the results are fairly insensitive 
both to the magnitude of the critical neck radius and to 
the bias potential strength m). and two parameters in a 
“suppression factor” that accounts for the decrease of the 
shell-plus-pairing correction with energy [3] . 

We have recently shown that the observed magnitude 
of the odd-even staggering can be directly correlated to 
the excitation energy above the outer part of the calcu¬ 
lated potential-energy surface [T^ ■ Therefore we suggested 
the BSM model could describe odd-even staggering if a 
potential-energy model were developed that accounts for 
how the individual nascent fragment properties are ex¬ 
pressed in the calculated potential-energy surface m- 

Here we propose a model for the potential-energy sur¬ 
face in which the properties of the individual fragments 
gradually emerge as the scission configuration is approached 
and specify the full details this proposed model. To treat 
odd-even staggering we add to the potential energy, as is 
customarily done since the dawn of nuclear mass calcu¬ 
lations [nun], a pairing contribution A to emerging odd 
fragments. However, since we start our trajectories at the 
ground-state shape of the even-even parent nuclei, where 
there is no odd pairing effect, only a fraction of the pair¬ 
ing delta of the fully formed final fragments is added in 
the initial stages of division. This fraction grows with de¬ 
creasing neck diameter, see below for further discussion 
and specification. 

We also introduce a correlated transfer of paired nu¬ 
cleon configurations. Such correlated nucleon transfers, 
which are different from sequential transfers, are often 
identified in nuclear reaction experiments, for a discussion 
see the presentation in Ref. m and the many references 
cited therein, for example Refs. [Min]- 

It turns out that we can incorporate our new ideas in 
calculations of charge yields, with minimal modifications 
of existing computer codes. Our first results on charge 
distributions are presented below. 


2 Shape parameterization and notations 

To describe the nuclear shape we use the three-quadratic- 
surface (3QS) parameterization. It was introduced almost 
50 years ago |18] . It is much more cumbersome to deal 
with than, say, a multipole expansion such as the /3 pa¬ 
rameterization m, but is used to allow a realistic descrip¬ 
tion of shapes of a fissioning nucleus up to and including 
division of the single shape into separated fragments. Par¬ 
ticularly noteworthy is that the emerging fragment shapes 
can be deformed spheroids or exact spheres. The latter is 
of special importance because it allows the extra binding 
associated with spherical doubly-magic nuclei to be accu¬ 
rately calculated. More details are found in Ref. (TS]; the 
discussions of Figs. 1 and 2 there are particularly infor¬ 
mative. How we design our potential-energy calculations 
in terms of this parameterization is detailed in Refs. [501 
nni; here we briefly summarize a few essential details. 



Fig. 1. Shape generated by the three-quadratic-surface pa¬ 
rameterization. Different colors distinguish between the shape 
sections generated by the three expressions in Eq. O- 

In the (3QS) parameterization the shape of the nu¬ 
clear surface is specified in terms of three smoothly joined 
portions of quadratic surfaces of revolution. They are com¬ 
pletely specified [TH] by 


Pizf 



< 02 ^ 






/i)^ , h - Cl < z < zi 
hY , < z <l 2 +C 2 

hf , Zi < Z < Z2 


( 1 ) 


Here the left-hand surface is denoted by the subscript 1, 
the right-hand one by 2 and the middle one by 3. Shapes 1 
and 2 are spheroids, for which c is the semi-symmetry axis 
length, a is the semi-transverse axis length, and I specifies 
the location of the center of the spheroid. The middle body 
may be a spheroid or a hyperboloid of one sheet, for which 
C 3 is imaginary. At the left and right intersections of the 
middle surface with the end surfaces the value of z is zi 
and Z 2 , respectively. Surfaces 1 and 2 are also referred to as 
end bodies and, alternatively, nascent fragments. A shape 
generated by the parameterization in Eq. o is shown in 
Fig.H 

In our calculations we use five shape parameters: elon¬ 
gation in terms of quadrupole moment Q 2 , left and right 
fragment spheroidal deformations efi and ef 2 , neck diam¬ 
eter d and mass asymmetry a^. These five parameters 
completely specify the shape for which the potential en¬ 
ergy is calculated and have been extensively discussed in 
Ref. [To]. Also, they completely exhausted the shape space 
available to this parameterization (except, trivially, for in¬ 
cluding larger values in each of the five shape parameters). 
For our studies here we need to revisit the asymmetry vari¬ 
able ttg which is directly connected to the asymmetry of 
nuclear shape (and the separated fragments): 

a-^ci U 2 O 2 

Ckg o o • 

afCi -I- 0202 
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This is equivalent to 


_ - Mi 

“ Mi + Mi 


(3) 


where and Mi are the volumes inside the end-body 
quadratic surfaces, were they completed to form closed- 
surface spheroids and “c” in the notation is to clarify 
we are referring to quantities of the compound single¬ 
shape system. To avoid the introduction of a large number 
of equivalent concepts we use somewhat interchangeably, 
“mass”, “volume” and nucleon number A. 

Our random-walk tracks end at a point where the neck 
radius of the nuclear shape is quite well developed, namely 
at 2.5 fm (compare to the radius of a spherical ^'^°Pu nu¬ 
cleus which is 7.2 fm), so the shape has almost reached 
the configuration of separated fragments. The neck radius 
is actually smaller than the radius of for which the ra¬ 
dius is 2.92 fm, and which is the lightest nucleus to which 
we have applied the macroscopic-microscopic method. Al¬ 
though Mi and Mi for shapes with well-developed necks 
can be expected to be close to the final fragment masses 
we cannot directly compare Mi and to the observed 
fission fragment masses Mf and M| (where the superscript 
“s” indicates we refer to the separated fragments) because 
the former do not quite sum up to the total nuclear vol¬ 
ume or mass A. However, by scaling Mi and Mi so that 
their sum adds up to the total mass number A, we can 
directly relate the “mass asymmetry” of the compound- 
nucleus shape to the observed heavy and light fragment 
masses. We obtain trivially 

1 + a 

Mi = r,Mi = 

Mi = r^Mi = A ^ and (4) 

A 

“ Mi + Mi 


where rg is a scaling factor for a nucleus with A nucle¬ 
ons. The scaling is equivalent to a assigning the mass in 
the neck region to the left and right nascent fragments 
in proportion to the respective volumes of these nascent 
fragments, were they completed to closed spheroids. The 
amount of matter involved is in the range of 10-20 nu¬ 
cleons. It is obvious how the same definitions apply to 
the proton and neutron numbers. The symbols Z and N 
are used for proton and neutron numbers of the fissioning 
compound system. We then have: 


Zi = Z-Zi and Ni = N- Ni (5) 


Therefore we usually use only Zf and Nf when we explic¬ 
itly refer to the number of fragment protons and neutrons, 
the other fragment proton and neutron numbers are then 
also specified. To calculate the yield as a function of both 
fragment proton and neutron number obviously requires 
that the previous 5D model (which always assumed that 
the proton to neutron ratios in both fragments were equal 
to the proton to neutron ratio in the compound system) 


be generalized to 6D so that the yield Y^Zf, Nf, Z^, Ni), 
which is a function of both proton and neutron asymme¬ 
try is obtained. We also realize that to describe odd-even 
effects requires that we space the grids we will use in terms 
of integer spacing of Zf and TV® in some fashion which we 
will now introduce. 


3 Potential energy versus shape and nascent 
fragment proton and neutron number 

The method to calculate the two-dimensional fission-fragment 
yield Y {Z, N) function that we now introduce is perhaps 
most transparently explained by occasionally referring to 
specific aspects of the computer code used to calculate, in 
the macroscopic-microscopic method, potential-energies as 
functions of shape. In our current potential-energy model 
and code we obtain a total potential energy for a specific 
compound nucleus (Z, N) and a specific “shape” as a sum 
of a macroscopic energy (given by a liquid-drop type ex¬ 
pression) and a microscopic shell-plus-pairing correction: 

Epot{Z, N, shape) = E^^dZ, N, shape) 

+A47(Z, (TV), shape) (6) 

+A4-‘((Z),iV, shape) 

where AI, shape) is the macroscopic energy and 

is the proton shell-plus-pairing correction and the 
final term the neutron shell-plus-pairing correction. In a 
calculation for a specific compound system and a specific 
shape these shell-plus-pairing corrections can trivially be 
individually tabulated. To obtain, say, the proton micro¬ 
scopic correction single-particle levels are calculated nu¬ 
merically in a folded-Yukawa single-particle potential with 
a functional form derived from the nuclear shape mm- 
From these levels the shell correction is obtained by use 
of the Strutinsky method [22112^ and the pairing correc¬ 
tion through, in our case, the Lipkin-Nogami method, as 
detailed in m- Thus, the proton microscopic correction 
is independent of neutron number, except that the poten¬ 
tial radius and depth depend on both proton and neutron 
number, therefore we have used the notation (N) and {Z) 
to indicate a weak dependence. But for small variations 
of neutron number the effect on the proton microscopic 
correction is small and can be neglected, an important 
observation that we will make use of later. To show this 
insensitivity we have calculated the proton and neutron 
shell-plus-pairing corrections for the ground-ground state 
shape of ^ 08^8162 for the single-particle fields appropriate 
for this nucleus. We find for the proton and neutron shell- 
plus-pairing corrections -3.7023 and -5.3715 MeV, respec¬ 
tively. When we do the calculation for single-particle fields 
appropriate for ii|Cni 62 we find that the proton and neu¬ 
tron shell-plus-pairing corrections are -3.6305 and -5.3353 
MeV respectively. So when we change the proton number 
by 4 units and implement the corresponding effect on the 
neutron single-particle field the neutron shell correction 
changes by only 0.0362 MeV. The change in the proton 
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shell-plus-pairing correction is larger because proton field 
changes more than the neutron field when we change the 
proton number. 

In the fission potential-energy code the shapes that can 
be studied are described in terms of the three quadratic 
surfaces of revolution: two end spheroidal sections and a 
middle region that near scission is a hyperboloid of revo¬ 
lution m as discussed in the previous section, see Eq. © 
and Fig. [TJ In the computer code the equations for these 
shapes need to be specified; five independent shape param¬ 
eters are required for the shape specification. Historically 
the shape parameters a and a were used m- But to more 
directly visualize the shape, we now characterize the shape 
in terms of five equivalent shape parameters: quadrupole 
moment Q 2 , related to the elongation of the shape, neck 
diameter d, left fragment spheroidal deformation efi, right 
fragment spheroidal deformation ef 2 , and mass asymme¬ 
try Q!g. The transformations from these parameters to the 
parameters of the quadratic functions that generate the 
shapes in the code are very lengthy and non-linear. They 
are described in Ref. [10] . For our discussion here we need 
to know that 


Ml - Mi 
~ Ml + Mi ~ Ml + Mi 


(7) 


Rather than discussing the asymmetry in terms of nu¬ 
cleon number A we can use the above concepts and discuss 
the charge asymmetry; we have earlier assumed that the 
proton to neutron ratio is the same in both fragments. We 
now develop an approach to treat different ratios, that is 
if the proton and neutron numbers are Zf and iV® in one 
of the fragments and Z — Zf and N — TV® in the other 
fragment we will treat 

ys y 7 7^ 7^ 

— A — A ^ = — (8) 

Nl^ N ^ N - Nl Ni ^ ’ 

It follows from previously that when we discuss the asym¬ 
metry in terms of proton number we can write 


Zl - Zi 
Zl + Zi 


(9) 


When we calculate the potential energy for the compound 
system that should correspond to specific (in our case inte¬ 
ger) separated-fragment charge numbers we use Eq. m to 
define the asymmetry ag of the corresponding compound- 
nucleus shape for which we calculate proton shell-plus- 
pairing corrections and macroscopic energies which are 
needed in our model. We need additional terms to describe 
how the macroscopic energy changes when we allow differ¬ 
ent proton to neutron ratios in the two fragments, mainly 
a symmetry-energy effect. We will discuss how to obtain 
this effect below. Correspondingly we can define the asym¬ 
metry ag for neutrons so that it corresponds to integer 
splits of neutron number and tabulate the calculated neu¬ 
tron shell-plus-pairing corrections. Below we specify how 
these results serve as the starting point to obtain the po¬ 
tential energy for ratios between the proton and neutron 
numbers that are different in the two fragments. 


We discussed above why the shell-plus-pairing correc¬ 
tions for protons and neutrons can be calculated indepen¬ 
dently of each other, to a very high degree of accuracy. 
Therefore we can write 


Epot{Z, N, Q2,d, Cfi, ££2, Zl, Nl) — 

Tfmac {Z,N,Q2,d,en,ef2,Zl,Nl) 

+ N, Q2,d, efi, ££2, agiZD) (10) 

+ N, Q2, d, ££1, ££2, ag(fVD) 

T Fodd 

Therefore, to obtain the total shell-plus-pairing correc¬ 
tions for any fragment split [Zf, Nl) we calculate and tab¬ 
ulate the proton shell-plus-pairing corrections for a grid 
in ag corresponding to integer Zf (and the corresponding 
Zi), obtained from Eq. We calculate the neutron shell- 
plus-pairing correction for a different spacing in ag corre¬ 
sponding to integer spacing in Nl. Thus the 6-dimensional 
shell-plus-pairing correction for any mass split (Zf,iVf) is 
the sum of two 5-dimensional functions. 

To obtain E^naclZ, N, Q 2 , d, ££i, ££ 2 , Zf, Nl) we proceed 
in several steps. First, when we run the code to calculate 
and tabulate the proton shell-plus-pairing correction for 
integer values of Zf we also tabulate the macroscopic en¬ 
ergy. It will then be obtained for non-integer values 

NI=N X ^ ( 11 ) 

Zi 

of Nl, because the asymmetry variable ag was chosen to 
correspond to integer values of Zf. 

Thus we have tabulated 

E^,,{Z, N, Q2,d, ££i, ££2, Zl, Nl) (12) 

where we need to remember that here Nl corresponds to a 
non-integer value because the asymmetry of the shape was 
chosen to yield integer Zf. The superscript “t” stands for 
“tabulated”. We need this tabulated value as one term in 
the macroscopic energy-model expression we now develop. 

But we need to obtain the macroscopic energy for (sev¬ 
eral different) integer iV®. It would be difficult to obtain 
such macroscopic energies by developing a model that in¬ 
tegrated across the compound nuclear shape for a config¬ 
uration with variable proton and neutron densities across 
the shape. But we now pose that we will get a sufficiently 
accurate model by considering changes in the macroscopic 
energy relative to the tabulated macroscopic energy we 
discussed in Eq. m- These changes are mainly due to 
changes in the symmetry energies, with considerably smaller 
contributions from other effects. We can obtain those by 
suitable consideration of changes in the macroscopic en¬ 
ergy of separated fragments. These we calculate as changes 
in sum of the macroscopic energy of separated spherical 
fragments. Therefore we calculate the sum of the spheri¬ 
cal macroscopic energies for two separated nuclei for the 
specific fixed Zl (related to Eq. (|T^ i and for a number of 
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Table 1. Spherical fragment macroscopic energies and their 
sums for various fragmentations of the compound system 
leading to fragment charges 52/40; only the neutron fragmen¬ 
tations vary, the columns labeled En and E {2 correspond the 
first and second term of line two in Eq. GSl). The lowest sum is 
obtained with the proton to neutron ratio as equal as possible 
in the two fragments and to that of the compound system. The 
line corresponding to this division is indicated by a “C” in the 
last column. 


z? 

Nf 

En 


m 

E{2 

Efi + Ef2 

52 

96 

-15.95 

40 

48 

-84.661 

-100.608 

52 

94 

-26.35 

40 

50 

-87.387 

-113.741 

52 

92 

-36.10 

40 

52 

-88.671 

-124.770 

52 

90 

-45.16 

40 

54 

-88.592 

-133.751 

52 

88 

-53.51 

40 

56 

-87.222 

-140.731 

52 

86 

-61.12 

40 

58 

-84.626 

-145.748 

52 

84 

-67.97 

40 

60 

-80.868 

-148.839 

52 

82 

-74.03 

40 

62 

-76.004 

-150.032C 

52 

80 

-79.26 

40 

64 

-70.088 

-149.348 

52 

78 

-83.64 

40 

66 

-63.170 

-146.807 

52 

76 

-87.12 

40 

68 

-55.298 

-142.422 

52 

74 

-89.68 

40 

70 

-46.515 

-136.199 

52 

72 

-91.28 

40 

72 

-36.862 

-128.143 

52 

70 

-91.87 

40 

74 

-26.377 

-118.250 

52 

68 

-91.42 

40 

76 

-15.097 

-106.515 

52 

66 

-89.87 

40 

78 

-3.055 

-92.926 


different N^: 

EZl{Zl,N,,Z-Zl,N-N,) = 

N,) + E^tiZ - Zl N - N,) (13) 

No odd-particle pairing effects should be included here; 
those are treated as discussed below. This function is tab¬ 
ulated for Z\ = 52 in Table 1 for fission of We 

note that in this integer-spaced grid the minimum energy 
occurs for a split where the sum of Z/N ratios in the 
two fragments is as close as possible to 2 x Z/N of the 
compound nucleus. The line corresponding to this split is 
indicated by a “C” at the very right. 

We pose that the macroscopic energy for any fragment 
division [Zf, Nf) in the fissioning system is given as 

EniEiC {Z,N,Q2,d,enNf2,Zl,N!) = 

^mac {Z,N,Q2,d,en,ei2,Zl,Nl) 

+ E^^^,{ZIN,,Z-ZIN-N,) 

- E^^iZl, Z - Zl, N - (14) 

where the last term is calculated by interpolation in the 
table corresponding to Eq. m- As a specific example 
we discuss a fragment division where the charge split is 
52/40. Then Nf = 81.39. We tabulate the sum in Eq. 
(fT51) as the right column in Table [T] and plot this sum 
for the specific charge division in our example in Fig. [51 


Heavy-Fragment Neutron Number 
90 88 86 84 82 80 78 76 74 



Light-Fragment Neutron Number 

Fig. 2. Sum of separated-fragment macroscopic energies. 

As an example that we can now, in our model, calcu¬ 
late the macroscopic energy for any (Zf, Nf). To illustrate 
this we continue to discuss our specific example. Equation 
m is a complete specification of the method. Suppose we 
want to calculate the macroscopic energy for the partic¬ 
ular case £'niac(92,144, Q 2 , d, Cfi, ef 2 , 52,64). Then we ob¬ 
tain the first term in the right member of Eq. m from our 
tabulated function in Eq. m- For our choice of Nf = 64 
the second term is given by the energy at the upper hori¬ 
zontal line and the third term by the energy at the lower 
horizontal line in Fig. [51 Thus in this example we obtain 
the macroscopic energy as a sum of the tabulated macro¬ 
scopic energy plus AE indicated in Fig. |2l 

One may argue that since some of the fissioning shapes 
involve deformed nascent fragments the terms in Eq. (US 
should be evaluated for the corresponding deformed shapes. 
But for the differences we consider here it only makes a 
minuscule difference. Let us choose ^°^Zr and ^°®Zr in Ta¬ 
ble |T] as an illustrative example. The energy difference for 
the spherical shapes AEsph tabulated is 

AE^pi, = -63.108 - (-75.973) = 12.865 (15) 

We have evaluated the macroscopic energies for deformed 
shapes with £2 = 0.5, which is the largest deformation for 
emerging fragments in our specified deformation grid, see 
Ref. [TU] and obtain 

AEdef = -53.262 - (-66.108) = 12.846 (16) 
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Thus, although the absolute energies change considerably 
the effect on their differences is completely negligible, only 
0.019 MeV. 

Finally we need to account for odd pairing effects in 
the emerging fragments. When the fragments are fully de¬ 
veloped with zero neck radius of the compound system a 
common assumption is that for odd-odd splits the extra 
odd contribution to the energy should be Eodd = 2 x A 
where A is the pairing-gap parameter, chosen as Z\ = 1.0 
MeV. With a non-zero neck radius the effect of pairing of 
the emerging odd splits would be smaller; for very com¬ 
pact shapes with no obvious neck there should be no effect. 
This leads to the following prescription for the odd energy 
of the fissioning system 

Eodd = 2x Ax {Bw - 1)'= odd Zf, 

£^odd = 0 even 

where, with the choice /c = 1, (i?w — 1)^ is the shape 
dependence of the Wigner term in our potential-energy 
model. The shape factor i?w is 1 for a shapes with no neck 
and increases continuously, as the neck develops, to 2 for 
separated fragments. It is necessary to postulate such a 
shape dependence because the macroscopic energy of two 
separated fragments contains two Wigner terms, the orig¬ 
inal system only one. Without such a shape dependence 
a discontinuity of the order of 10 MeV would occur at 
scission of actinide nuclei. That is, if we calculated the 
energy for a deforming nucleus up to the scission point 
we would at scission obtain a 10 MeV lower energy than 
if we calculated the energy of two approaching separated 
fragments. A pedagogical figure illustrating this and the 
necessity of this shape dependence is in Ref. E], Fig. 1. 
Since we need a realistic potential-energy surface in the 
scission region we do need to consider these issues (which 
have in many investigations been ignored). The compre¬ 
hensive discussion of the shape-dependence of the Wigner 
term in Ref. [3^ carries directly over to how the effect of 
the pairing A increases as the neck becomes more narrow. 
There is no known derivation of the Wigner shape depen¬ 
dence so it is just postulated, but with consideration of 
its limiting behavior |25j . The power constant k, which 
we introduced here to allow some sensitivity studies, gov¬ 
erns how early in the division process the character of the 
two fragments causes the “second” Wigner term, or in our 
case, the odd pairing effect, to manifest itself. Below we 
will present sensitivity studies on the shape dependence 
and on the magnitude of the pairing delta. 


4 Application to Charge Asymmetries 

A calculation of the complete (2D) Y{Z,N) yield distri¬ 
bution based on the above model would lead to much 
more complex calculations compared to our current 5D 
implementation Dills], because of its 6D nature and as¬ 
sociated vastly increased storage requirements. But, as a 
test of the above approach, we can study many of its as¬ 
pects by looking at the odd-even staggering in charge dis¬ 
tributions without treating the neutrons separately. We 
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Fig. 3. Two-dimensional slice in the 5D space schematically 
illustrating possible candidate points for the next trajectory 
step. The smaller square with gray-filled circles inside, indi¬ 
cates the limit originally chosen for next step candidate points. 
The circles partially filled with gray correspond to transfer of 
paired configurations. 


obviously then have to calculate the potential energy for 
field asymmetries Og that correspond to integer spacing 
of the proton number Zf. For for each combination 

of the other 4 shape variables [10], we calculate the po¬ 
tential energies for asymmetries corresponding to charge 
splits 47/47, 46/48, 45/49 ..., and accordingly for other 
elements. This corresponds to “averaging” or “summing” 
over the neutron variable, a procedure we assume has lim¬ 
ited effects on the charge distributions obtained relative to 
summing over N a calculated complete 2D Y{Z, N) dis¬ 
tribution. Thus, we obtain in strict analogy with our pre¬ 
vious results, on a discreet grid the 5D potential-energy 
matrix Epotih, I 2 , h, h, h) where Ii corresponds to the 
quadrupole moment (or elongation), I 2 to the neck radius, 
I 3 to spheroidal deformations of one of the emerging frag¬ 
ments, I 4 to spheroidal deformations of the other emerg¬ 
ing fragment, 15 to charge asymmetries ZIIZ 2 , where the 
charge numbers are integers. However, since the calculated 
energy, ifcomp, does not contain contributions from pair¬ 
ing effects in the emerging fragments we add the shape- 
dependent odd enhancement according to Eq. (IT71) . It turns 
out very few modifications of the random walk code are 
required for this calculation. 

In the BSM model we find the yield distributions by 
generating paths through the potential-energy matrix as 
follows. We usually start the path at the ground state 
(black dot in the schematic 2D Fig. [3]). We then select 
randomly one of the 242 surrounding points (circles with 
gray interior in Fig. D only 8 of them in this schematic 
representation) as a candidate for the next point on the 
track. Suppose the energy difference between this point 
and the current point is AV. Then, if AV < 0 the can¬ 
didate point becomes the next point (and current point 
“at the next throw of the dice”) on the path. However, 
if AV > 0 this outcome is only selected with probabil¬ 
ity P = exp{—AV/T) where T is the temperature; full 
details are in Refs. m- When the critical neck radius 
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Fig. 4. Calculated and measured charge yields for neutron- 
induced fission of The top panel represents the original 
model, the middle panel has 2xZi added to the potential energy 
in the matrix locations corresponding to odd charge splits. In 
the bottom panel steps corresponding to two-proton changes 
in asymmetry are allowed. 

Co = 2.5 fm is reached the walk is terminated and the 
asymmetry of the shape recorded. Each yield curve is 
based on 20000 tracks. We present results for four re¬ 
actions: thermal neutron-induced fission of 

and photon-induced fission with energies centered 
around 11 MeV for in Figs. SI O 0 and [S] respec¬ 
tively. The experimental data for the (n,f) reactions are 



30 40 50 60 

Fragment Charge Number Z, 


Fig. 5. As Fig. SI but the experimental data are for ( 7 ,f) reac¬ 
tions leading to E* « 8 — 14 MeV; they include contamination 
from fission after ln(« 15%) and 2n(Ri 5%) emissionS]. The 
calculations are for fission of at E* = 11 MeV. 

from Ref. [55], the ( 7 ,f) data from Ref. S]. The top frame 
in each of the four figures is with the original model with 
no pairing effect added. The middle frames are based on 
BSM in the 5D potential modified according to Eq. (SZD. 
We find little staggering in the calculated curves, although 
we implemented the “standard” explanation for odd-even 
staggering: for odd-odd splits the potential is on average 
2 X A (for zero neck radius; less in earlier stages of the 
fission process) higher than the potential for even splits. 
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Fragment Charge Number Z, 


Fig. 6. As Fig.Slbut for 

But it is clear that the original formulation of the BSM 
model would never be able to describe odd-even stagger¬ 
ing, even after the addition of a pairing effect to the cal¬ 
culated potential energy. To illustrate why let us look at 
Fig. 1 and specifically at .Z = 38 on the experimental 
curve. Let us assume that Z = 38 is the asymmetry of 
the current point on our evolving track. Because we only 
consider next-neighbor gridpoints as candidate points for 
the next point on the path, we have to populate Z = 39 
which has a very low yield, to subsequently populate the 
high-yield Z = 40 point, for example. To obtain a pro¬ 
nounced staggering with these model features is impossi¬ 



30 40 50 60 

Fragment Charge Number Z, 


Fig. 7. As Fig.Slbut for ^^Pu. 

ble. But as the shape evolves in the asymmetry direction 
and two levels cross it is reasonable that an alternative 
to breaking a pair and transferring only one proton be¬ 
tween the evolving fragments is that a paired two-proton 
configuration could be transferred in one step. We have 
implemented this possibility in the BSM model by also 
allowing Z — 2 and Z -|- 2 as next track-point candidates 
(shown as circles partially filled with gray in Fig. [3]). As al¬ 
luded to above transfer reactions indicate that correlated 
transfer of nucleon pairs are common see Refs. [T51IT61IT7] . 
In our current treatment transfer of either a paired con¬ 
figuration or breaking of a pair and transfer of one proton 
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30 40 50 60 

Fragment Charge Number Z, 


Fig. 8. Calculated charge yield for neutron-induced fission of 
^^°Pu in the BSM model for different assumptions about pair¬ 
ing and next-step grid-point candidates. In the top panel no 
pairing effect is added to the potential energy for odd splits, but 
one or two steps in asymmetry is implemented. In the middle 
panel we add a pairing effect to odd splits and permit as next 
step candidates one or two grid points in asymmetry and both 
of the fragment deformations. In the bottom panel we have 
also allowed one or two steps in the elongation coordinate Q 2 . 


can occur. Typical excitation energies at the end of the 
asymmetric tracks, see Ref. [12] are 7.7-11.6 MeV, and 
lower earlier in the shape evolution. Our consideration of 



30 40 50 60 

Fragment Charge Number Z, 


Fig. 9. Study of yield sensitivity to the magnitude of the pair¬ 
ing A and the onset of fragment pairing effects. Except for the 
very small pairing effects on which the bottom panel is based 
the results are quite robust and similar to the standard result 
in the bottom panel in Fig. [T] 

paired configurations is quite consistent with the results 
of Ref. [27], where in one example at 8.4 MeV excitation 
the proportion of paired configurations is 36%. 

The calculated yields with transfers of correlated pairs 
allowed as next track-point candidates are in the bottom 
panels in Figs. HHZl Staggering is now obvious in the cal¬ 
culated results and in close agreement with the experi¬ 
mental data. In the calculated results in Figs. |3H7]it may 
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seem that the crucial generalization that we introduced 
to describe the odd-even staggering is not the first step 
we took, namely the addition of 2 x Z\ to the odd charge 
splits, but the second step in which we permitted a change 
in Z of two units, corresponding to a transfer of a paired 
proton configuration. But what if we had implemented 
this as a first step? We have investigated this possibility 
and show in the top panel of Fig. [5] the calculated charge 
yield for neutron-induced fission of ^'^^Pu with no pairing 
energy added to the odd charge splits, but with both one 
and two steps in Z permitted as next candidate points on 
the random-walk trajectory. There is no odd-even stag¬ 
gering in the calculated curve which is extremely similar 
to the calculated yield in the top panel in Fig. [T] It is in¬ 
teresting to observe that although allowing both one and 
two steps in Z effectively corresponds to increasing the 
speed of motion in asymmetry, or equivalently increasing 
the distance between grid points, there is little change in 
the calculated yield curve. In the middle panel of Fig. |5] 
we have added (fractions of) 2Z\ to the potential energies 
for odd charge splits, allow one or two steps in asymmetry, 
but also allow one or two steps in the fragment deforma¬ 
tion shape variables. There is little difference compared to 
the bottom panel in Fig. [T] This again shows, as was ear¬ 
lier pointed out [5] that the calculated yield curves in the 
BSM model do not depend sensitively on most changes in 
the deformation grid. In the bottom panel of Fig. |8] we 
have furthermore allowed one or two steps in the elonga¬ 
tion variable Q 2 - In this case there is a noticeable effect 
on the calculated yield. Now the calculated distribution 
is slightly wider than the experimental results. In Ref. [5] 
we showed that tripling the number of points in Q 2 led to 
a calculated distribution that was narrower than the ex¬ 
perimentally observed one. Obviously, extreme changes in 
the grid will influence the calculated results. For example 
if we were to use only three grid points in the elongation 
variable we would not obtain realistic yields. 

Finally we study, for ^^°Pu, the sensitivity to variations 
of the postulated shape dependence that governs the on¬ 
set of odd-even effects on the potential energy and to the 
magnitude of the pairing A. The results are shown in Fig. 
E The top panel shows the effect of increasing the pairing 
A by 20% relative to our standard assumption of = 1.0 
MeV. The calculated curve is little different for the result 
in the bottom panel of Fig. [T] Therefore the exact choice 
of Z\ is a non-issue in this first study of odd-even stagger¬ 
ing. In larger systematic studies a well-defined prescription 
should obviously be introduced. The middle and bottom 
panels show the results for different forms of the shape 
factor (i?w ~ 1)^ which governs the rate with which the 
effect of the odd-even pairing A in the fragments affects 
the calculated potential energy as scission is approached. 
The quantity is close to 1.5 at our selected scission 
configuration. Therefore, in the middle panel the factor 
(i?w — 1)°'® is 0.7, so that for the odd configurations 1.4 
MeV is added, less earlier in the division process. The stag¬ 
gering here is only very slightly larger than in our standard 
calculation in the bottom panel of Fig. [7] where 1.0 MeV is 
added in the scission region. In the bottom panel of Fig. IH] 


the shape factor (Rw — 1)^ comes out to 0.25 so that only 
0.5 MeV or less is added to the odd-odd divisions. Here 
the staggering is much reduced, as might be expected. 
However from these sensitivity studies we conclude that 
the results are quite robust for reasonable variations of 
assumptions about the onset of fragment character on the 
potential energy as well as to the magnitude of the pairing 
Zi. 

In summary we have shown that to describe odd-even 
staggering in the BSM model we must add simultane¬ 
ously two effects: (1) odd-even effects on the calculated 
potential-energy surface and (2) allow transfers of corre¬ 
lated paired proton configurations; they work together in 
the development of odd-even staggering. 

As future “perspectives for the next decade” we antic¬ 
ipate that to develop more accurate descriptions we need 
to 

— implement the extensions we discuss here into a com¬ 
puter model framework and calculate the full Y{Z, N) 
fission-fragment yield distributions (a two-dimensional 
function of neutron and proton number), 

— treat from more basic principles how the number of 
paired configurations decrease with excitation energy 
(see Ref. m for a discussion) which should influence 
the probability with which a point corresponding to 
two-nucleon transfer is chosen as the next candidate 
point on the trajectory, 

— calculate the damping of shell effects based on actual 
single-particle structure rather than use a parameter¬ 
ized approach. 

— understand issues related to the deformation grid. Clearly 
the calculated yields do depend on the selection of the 
grid. To take an obvious example, were we to have only 

3 grid points in the elongation direction we would not 
obtain any realistic yields. However, in the extensive 
sensitivity studies in Ref. [5] it was shown that the 
yields were remarkably insensitive to many types of 
grid changes, as is also shown by the results of Figs. [5] 
and ini above. 
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